Evolution of Quantum Discord and its Stability in Two-Qubit NMR Systems 
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We investigate evolution of quantum correlations in ensembles of two-qubit nuclear spin systems 
via nuclear magnetic resonance techniques. We use discord as a measure of quantum correlations 
and the Werner state as an explicit example. We first introduce different ways of measuring discord 
and geometric discord in two-qubit systems and then describe the following experimental studies: 
(a) We quantitatively measure discord for Werner-like states prepared using an entangling pulse 
sequence. An initial thermal state with zero discord is gradually and periodically transformed into 
a mixed state with maximum discord. The experimental and simulated behavior of rise and fall 
of discord agree fairly well, (b) We examine the efficiency of dynamical decoupling sequences in 
preserving quantum correlations. In our experimental setup, the dynamical decoupling sequences 
preserved the traceless parts of the density matrices at high fidelity. But they could not maintain 
the purity of the quantum states and so were unable to keep the discord from decaying, (c) We 
observe the evolution of discord for a singlet-triplet mixed state during a radio-frequency spin-lock. 
A simple relaxation model describes the evolution of discord, and the accompanying evolution of 
fidelity of the long-lived singlet state, reasonably well. 
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I. INTRODUCTION 



Since its introduction by Schrodinger, entanglement 
has remained an extensively studied and yet a mysterious 
aspect of quantum theory. Entanglement appears as a 
by-product of the quantum formalism that assigns prob- 
ability amplitudes to physical states and lets them exist 
in coherent superpositions. Although, it runs counter 
to the human intuition gained through experiences with 
classical systems, experimental evidence has consistently 
favored the existence of such superposed quantum states. 
Entanglement was thought to be an indispensable re- 
source for quantum information processing, which can 
outperform the corresponding classical information pro- 
cessing. Indeed, various quantum algorithms that exploit 
entanglement have been proposed and successfully tested 
[H . Separable (i.e. not entangled) states were considered 
insufficient to implement quantum information process- 
ing. That belief has changed since Ollivier and Zurek 
as well as Henderson and Vedral [S] independently intro- 
duced a new measure of non-classical correlations named 
^discord\ Discord is based on the measure of mutual 
information between two parts of a system. It can be 
put in one-to-one correspondence with entanglement for 
pure states, but, unlike entanglement, it can be nonzero 
for separable mixed states [3] . 

There have been long debates on the necessity of en- 
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tanglement for quantum information processing. For ex- 
ample, the purity of typical spin systems used in nuclear 
magnetic resonance (NMR) experiments is too small to 
exhibit entanglement; nevertheless, NMR quantum in- 
formation processing is considered an efficient test bed 
for quantum algorithms [^-l6|. Furthermore, Knill and 
Laflamme proposed an algorithm, called DQCl, which 
estimates the trace of any unitary matrix faster than 
any known classical algorithm ^]. Datta et al. showed 
that entanglement in the DQCl algorithm is vanishingly 
small, and it further decays with increase in the number 
of qubits Q. They also showed that the DQCl algorithm 
involves nonzero discord. Thus, our present notion of 
quantum speed up is tied to discord rather than entan- 
glement [9j. For a more detailed review on discord, see 
Ref. [1^. 

Discord happens to be one of the many quantities that 
can measure the nonclassicality of a given quantum sys- 
tem. Its standard definition equates discord to the dif- 
ference between two classically equivalent forms of mu- 
tual information. Due to the difficulty in measuring this 
difference, Dakic et al. proposed an alternate quantity 
called ^geometric discord\ It is the distance between the 
given quantum state and the closest classical state and is 
easier to quantify than discord [ll|. We look at both of 
these quantities in our study. 

Several experiments to demonstrate quantum correla- 
tions in liquid-state NMR systems have been performed 
in recent years, e.g., by measuring a suitable witness op- 
erator [I4I, by measuring discor d llSll . and by evaluat- 
ing the Leggett-Garg inequality [l^. Measurements of 
discord in mixed states have also been performed using 



optical systems [15| and quadrupolar NMR systems [16| . 
Here we report an experimental study of time evolution 
of discord in NMR systems. After preparing a two-qubit 
Werner state, we study the accumulation of discord, ef- 
fects of dynamical decoupling sequences on it, and its 
decay due to decoherence. In Sec. II, we revisit some 
theoretical aspects of discord and describe different ways 
of measuring it for the Werner state. Then, in Sec. Ill, 
we present the experimental details and discuss the re- 
sults. We conclude in Sec. IV with some inferences from 
our analysis. 



II. THEORY 



Discord 



system, and is given by 



I{A : B) = n{A) + H{B) - H{A, B) . 



(5) 



This expression can be intuitively understood with the 
help of Fig.[TJ On the right-hand side, the first two terms 
quantify the information content of subsystems A and B, 
respectively. So the information common to both the 
subsystems is counted twice. Subtracting the informa- 
tion content of the combined system then gives the com- 
mon (or mutual) information. The result is clearly sym- 
metric, i.e., X{A : B) = X{B -.A). A classically equiv- 
alent form of mutual information, also shown in Fig. [TJ 
is 



J{A:B) = n{B) - n{B\A) 

- H{B)-Y,PfnB\a 







(6) 



Conditional Entropy: In classical information theory 
the amount of information contained in a random vari- 
able X is quantified as the Shannon entropy. 



"HCX) = -^p^l0g2Px 



(1) 



where p^. is the probability of occurrence of event X. 
When T-L{X) ~ 0, the random variable X is completely 
determined and no new information is gained by mea- 
suring it. Hence, Shannon entropy can be interpreted as 
either the uncertainty before measuring a random vari- 
able or the information gained on measuring it. 

Consider a bipartite system containing two subsystems 
(or random variables), A and B. Conditional entropy of 
B quantifies the uncertainty in measurement of B when A 
is known, and is represented by 'H{B\A). Using classical 
probability theory, it can be expressed as 



'H{B\A) ^H{A,B)~H{A) , 



(2) 



where 'H(A, B) is the information content of the full sys- 
tem and HiA) is the information content of the subsys- 
tem A. An equivalent way of defining the conditional en- 
tropy is 



n{B\A)=J2p'^nB\a = i) , 



(3) 



which removes from the information content of subsystem 
B the conditional contribution that is not contained in 
subsystem A. 

In quantum information theory, the von Neumann en- 
tropy gives the information content of a density matrix. 



H{p) =-^Xx log2 A^ 



(7) 



where A^; are the eigenvalues of the density matrix p. 
Although the two expressions of mutual information, 
Eqs. ([5]) and ^, are equivalent in classical information 
theory, they are not so in quantum information theory. 
The reason for the difference is that the expression for 
mutual information given by Eq. ([SJ involves measure- 
ment and depends on its outcomes. Measurements in 
quantum theory are basis dependent and also change 
the state of the system. Henderson and Vedral j3| have 
proved that the total classical correlation can be obtained 
as the largest value of 



J{A : B) = H{B) - H{B\A) 

= H{B)-Y,P^H{B\a = i) 



!H{K) 5f(B) 



(8) 



where 



n{B\a ^ i) ^ -Y,p{h,\a,)\og2p{h,\a,) , (4) 



i5f(A|B) J(A:B) 5f(B|A) 



and p(bj\ai) is the conditional probability of occurrence 
of event bj given that event a^ has occurred. Unlike the 
definition in Eq. ©, the definition in Eq. ^ involves 
measurement of one subsystem of a bipartite system. 

Mutual Information: It is the amount of information 
that is common to both the subsystems of a bipartite 



5^A,B) 

FIG. 1: The Venn diagram representing total information 
'H{A,B), individual informations {H{A), H{B)), conditional 
information ('^(^1-8), 'H{B\A)), and mutual information 
T{A : B) = J{A : B) in classical information theory. 



where the maximization is performed over all possible 
orthonormal measurement bases illfl for A, satisfying 
'E.Uf = 1 and n^n^ = 5„n^ [13. Therefore, the non- 
classical correlations can be quantified as the difference 



D{B\A) = I{A 



B)-maxJ(A:B) . (9) 

{n°} 



Ollivier and Zurek named this difference 'discord^ [21- 
Zero-discord states or "classical" states are the states in 
which the maximal information about a subsystem can 
be obtained without disturbing its correlations with the 
rest of the system. 

Discord is not a symmetric function in general, i.e. 
D{B\A) and D{A\B) can differ. Datta ')M, has proved 
that a state pab satisfies D{B\A) = if and only if 
there exists a complete set of orthonormal measurement 
operators on A such that 



PAi3 = Ep^^^ 



PB\a= 



(10) 



When one part of a general bipartite system is mea- 
sured, the resulting density matrix is of the form given by 
Eq. (fTO| . Since the state rendered on measurement is a 
classical state, one can extract classical correlations from 
it. Thus, for any quantum state and every orthonormal 
measurement basis, there exists a classically correlated 
state. Maximization of J{A : B) captures the maximum 
classical correlation that can be extracted from the sys- 
tem, and whatever extra correlation that may remain is 
the quantum correlation. 



B. Evaluation of Discord 

Given a density matrix pab , it is easy to construct the 
reduced density matrices pA and Pb and then obtain the 
total correlation I{A : B) using the quantum analog of 
Eq.©. Maximization of J {A : B) to evaluate discord is 
nontrivial, however. The brute force method is to max- 
imize J{A : B) over as many orthonormal measurement 
bases as possible, taking into account all constraints and 
symmetries. For a general quantum state, a closed ana- 
lytic formula for discord does not exist, but for certain 
special class of states analytical results are available [1^ . 
For example, Chen et al. have described analytical eval- 
uation of discord for two qubit X states under specific 
circumstances [201. Luo has given an analytical formula 
for discord of the Bell-diagonal states that form a subset 
of the X states [21j. In our work, we evaluate discord 
using both the brute force method and the Luo method. 
We use Bell-diagonal states in our experiments, but our 
experimental preparation of the states is not perfect. The 
difference between the discord values obtained by the two 
methods then provides an estimate of the experimental 
imperfections. 

Extensive measurement m,ethod: This method involves 
measurements over extensive sets of orthonormal basis 



vectors and maximization of J{A : B). For measure- 
ment of a single qubit in a two-qubit system, we use the 
orthonormal basis 

{\u) = cos 6*10) +e*"* sin 6l|l) , \v) = sin6l|0)-e*'^cos6l|l)} , 

(11) 
and let cos 9 e [—1,1] and (p G [0, 27r) vary in small steps. 
For every choice of 9 and 0, we project the experimen- 
tal density matrix obtained by tomography along the or- 
thonormal basis. The postprojection density matrix is 



Z^ 



n'^puf = 



E 



pfU^<»PB\, 



(12) 



with p^ = Tr[Ilfp]. Discord is then obtained from the 
conditional density matrix pB\a=i as per Eqs. ([5]) and 

Strictly speaking, this method gives a lower bound on 
J {A : B), since the direction maximizing J{A : B) may 
not exactly match any of the points on the discrete (^, (j)) 
grid. Also, when the desired state is isotropic, e.g., the 
Werner state, the angular variation of J{A : B) provides 
an estimate of the inaccuracy in the state preparation, 
e.g., due to inhomogeneities and pulse imperfections. 

Analytical method for the Bell-diagonal states: As the 
name suggests, the Bell-diagonal states are diagonal in 
the Bell basis, given by 



1 



|V'±) = ^(|oi)±|io)) 



|0±) = -i=(|OO)±|ll)). (13) 



The generic structure of a Bell-diagonal state is pbd = 

Aiiv-)(V'-| + A2ir)(n + A3i0+)(0+i + A4iv+)(V'+|. 

With only local unitary operations (so as not to alter the 
correlations) , all Bell-diagonal states can be transformed 
to the form 



1 



Pbd 



E^ 



(14) 



where the real numbers rj are constrained such that all 
eigenvalues of pbd remain in [0, 1]. The symmetric form 
of Pbd also implies that it has symmetric discord, i.e., 
Dbd{B\A)^Dbd{A\B). 

Luo chose the set of measurement bases as {Fnjjyl'}, 
where IIJJ = \k){k\ are the projection operators for the 
standard basis states {k = 0,1), and V is an arbitrary 
SU{2) rotation matrix. A projective measurement yields 
the probabilities Po = Pi — 2 and an analytical formula 
for the classical correlation, 

max J(A : B) = {^^^ log2(l-^)+(^) ^og^il+r) 

(15) 
with r = max{|ri|, |r2|, Ira]}. 

For the Bell-diagonal states, the reduced density ma- 



trices are pA ^ Pb = 1/2, and the total correlation is 



quantum state and the nearest classical state, 



I{A: B) = 2 + Y^ A,: logs -^i 

i=l 

where the eigenvalues Xi of psD are 

Ai = (1-ri -r2-r3)/4 

As = (l-ri+r2 + r3)/4 

A3 = (1 + ri - 7-2 + r3)/4 

A4 = (1 + ri + 7-2 - r3)/4. 



(16) 



(17) 



Thus the analytical formula for discord is, as per Eq. ©, 

4 



Dbd{B\A) = 2 + ^Aaog2A, - 



1=1 
1 + r 



1-7- 



log2(l + r) . 



l0g2(l-7-) 



(18) 



For a Werner state of the form 
1-e 



Pwie) 



e\i;-){r 



(19) 



-e and r ~ e. The discord is then given by 

1 (l-e)(l + 3£) 

^^(^) - 4^"^^ (1 + e)^ 



lOE 



(l + 3e)3 



i2 ■ 






(20) 



This expression is plotted versus the purity e in Fig. [21 
together with the corresponding correlations I {A : B) 
and J {A : B). 

In practice, the experimental density matrix obtained 
by tomography is not necessarily Bell diagonal. We ob- 
tain I {A : B) as before, using Eq. ([5]). To extract 
the maximum value of J [A : B), we drop the the off- 
diagonal terms, keeping only the terms in Eq. (|14p. and 
use Eq. (jisp . In this procedure, discord is overestimated, 
whenever the actual direction maximizing J{A : B) is 
not in the Bell-diagonal state subspace. 



DG{B\A)^uAn\\p~x\ 



(21) 



where VIq represents the set of classical states, and 
||X-r||2 == Tr(X-r)2 is the Hilbert-Schmidt quadratic 
norm. Obviously, DG{B\A) is invariant under local uni- 
tary transformations. Analytical formula for computing 
geometric discord for an arbitrary Amxm ® Bnxn state 
of a bipartite quantum system is available [22] . Recently 
discovered ways to calculate lower bounds on discord 
for such general states do not require tomography and, 
hence, are friendlier experimentally |23ll24|. 

We follow the formalism of Dakic et al. [il| to obtain 
geometric discord for two-qubit states. The two-qubit 
density matrix in the Bloch representation is 

-.3 3 3 

i— 1 z— 1 'f'd—^ 

(22) 
where Xi and yi represent the Bloch vectors for the two 
qubits, and Tij = Tr[(p(tTi ® '^j))] are the components of 
the correlation matrix. The geometric discord for such a 
state is 

DG{B\A) ^ ^{\\x\\' + \\Tf - 7^^,,^) , (23) 

where HTp = Tr[T^r], and rymax is the largest eigenvalue 
of xx^ + TT^. 

For the Werner state, as already mentioned, Xi = Q = 
yi and T is a diagonal matrix with Ta — — e. Then 
||T|p = 3e^ and all eigenvalues of TT^ are e^, yielding 



DGw{e)^\{3e'~e')^'-. 



(24) 



This expression is also plotted versus the purity e in 
Fig. [2j Comparison with Eq. ((20)) reveals that dis- 
cord and geometric discord are proportional for low- 
purity Werner states. Also, the numerical difference be- 
tween Dwi^) and 2DGw{^) does not exceed 0.027 for all 
EG [0,1]. 

The Bloch parameters Xi, yi, and Tij provide a com- 
plete description of any two-qubit state. So tomographic 
measurement of these parameters determines the geomet- 
ric discord exactly by Eq. 



C. Geometric Discord 

Since the maximization of J{A : B) involved in calcu- 
lating discord is a hard problem, Dakic et al. introduced 
a more easily computable form of discord based on a ge- 
ometric measure [11]. For every quantum state there is 
a set of postmeasurement classical states, and the ge- 
ometric discord is defined as the distance between the 



III. EXPERIMENT 

We now describe experimental evaluation of discord for 
two-qubit NMR systems. We measured quantum correla- 
tions for two different samples, each forming a two-qubit 
system, under different circumstances. 
Sample 1 is ^H and ^'^C spins of ^■^C-chloroform [see 
Fig. [3^] dissolved in deuterated chloroform (CDCI3). 
Both ^H and ^'^C spins were on-resonant and the scalar 



coupling (J) between the two spins is 219 Hz. For the 
proton, the Ti and T2 relaxation time constants are 14.5 
and 5.7 s, respectively. For carbon, they are 21 and 0.25 s, 
respectively. 

Sample 2 is ^H nuclear pairs of 5-chlorothiophene- 
2-carbonitrile (see Fig. [6)d) dissolved in deuterated 
dymethylsulfoxside (DMSO-Dg). The chemical shift dif- 
ference (Ai/) and scalar coupling (J) between the two 
spins are 270 and 4.11 Hz, respectively. For each proton, 
the Ti and T2 relaxation time constants are 6.3 and 2.3 s, 
respectively. 

All the experiments were carried out in a Bruker 500- 
MHz NMR spectrometer at an ambient temperature of 
300 K. Precise radio-frequency (rf) gates for the exper- 
iments were synthesized by numerical optimizations as 
described in Refs. [2^, [2^. Decoherence in these systems 
is mainly due to the fluctuations in the local magnetic 
field, caused by random molecular motions in the pres- 
ence of intra- and intermolecular dipolar interactions and 
chemical shift anisotropics. Traces of paramagnetic im- 
purities also contribute to spin- relaxation 27 1. 

The Hamiltonian for a two-qubit NMR system, with 
spins I^ and I^ ^ can be written as 






J ■ 



(25) 



Here J^z = —fi {^aI^ + ^bIz) i^ the Zeeman Hamil- 
tonian, characterized by the Larmor frequencies to a and 
ujb , and J^j = 2ttJI^ ■ I^ is the indirect spin-spin cou- 
pling Hamiltonian. 

In thermal equilibrium at room temperature, kT is 
much larger than the Zeeman energy splittings. So the 
density matrix of a two-qubit system can be expanded as 

M 
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FIG. 3: (a) The molecular structure of ^■^C-chloroform, and 
(b) the pulse sequence for discord preparation. In (b), PFG 
is the pulse field gradient operation which destroys the coher- 
ences and retains the diagonal elements of the density matrix. 



The identity 1 represents a background of uniformly pop- 
ulated levels, and the traceless part Ji = S,{lf + ^I^) 
is known as the deviation density matrix. Only the trace- 
less part p is manipulated by unitary transformations in 
all NMR experiments. For protons in a magnetic field of 
strength 11.7 T at room temperature, the small dimen- 
sionless number ^ = hujA/kT sa 8 x 10~^. The discord 
for this size of purity is shown in the inset of Fig. [2] by a 
shaded area. 



A. Preparation of non-zero discord states 
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FIG. 2: Various correlations as functions of the purity e for 
Werner states of the form given in Eq. (|19|) . The inset shows 
the range of discord for the purity available in our NMR setup. 



We prepared the maximum discord state starting from 
the zero discord thermal equilibrium state (peq) as fol- 
lows. The pulse sequence in Fig. [3] was used to prepare 
discord between ^H and ^'^C spins of ^'^C-chloroform. For 
this system, ojc/^h ~ 1/4. An initial |00) pseudopure 
state was prepared using the spatial averaging method 
[23], as shown in the first part of Fig.[3ljb). The transfor- 
mations of the traceless Poq/^ under the spatial averaging 
pulse sequence are (as per standard notation _27] ) 

oA 



115] 

/^cos(15°)-/^sin(15°) 



4-^z 



// cos(15°) + 24^/f sin(15°) + i/f 
iTBtytGz 
I^ cos 15° cos 75° + 2I^lf sin 15° sin 75° 4 



IjB 



(27) 



This pseudopure state was converted into a Werner state, 
using the second part of the pulse sequence in Fig. ^h) 



with the delay e/(-nJ) = 1/(2 J): 



Ppp 



(4^ 



-'a; ' ^^ — X 



/f + 24^/f 

oB 

— X 

B otAtB 



i 90°^ 

+ 2,7 



oB 



PW 



; 180°-^, 90^ 
|(-24^/f-24^/f-2/,^/,^) 






(28) 



Comparing with Eq. ([T9| . we can see that the relevant 
purity parameter in this case is e = ^/8. 

We also varied 6 from to 27r in 13 steps and, for 
each delay, carried out toinography to measure the ex- 
perimental density matrix [29| . The corresponding simu- 
lated density matrices are obtained by assuming perfect 
pulses and carrying out a calculation similar to Eq. ([28| . 
It can be easily seen that for 6 values that are odd multi- 
ples of 7r/2, one obtains Bell-diagonal states. The fidelity 
i^ of a test density matrix ptcst relative to the Werner 
state is defined as [2a 



F = 



Ti-[Ptost • Pw] 



viMpL 



n-pw] 



(29) 



Fidelities of the experimental and the simulated density 
matrices, as functions of 9, are shown in Fig. Ufa). The 



discord for each value of 9 is obtained using the extensive 
measurement method as described in Sec. II-B. Both ex- 
perimental and simulated values of the discord are plot- 
ted in Fig. |li;b), in units of e2/ln2. The state at 6* = 
is related to the pseudopure |00) state by local unitary 
transformations and, therefore, has zero discord. Other- 
wise, for 9^0, nonlocal spin-spin interactions give rise to 
discord. For 9 equal to odd multiples of 7r/2, one obtains 
Bell-diagonal states with maximum discord. For 9 equal 
to TT, the delay equals the period of the scalar coupling, 
implying no transformation, and the discord is periodic 
thereafter. 

The nonzero discord in the above experiments indi- 
cate the intrinsic quantumness of nuclear spin systems 
even at high temperatures. Among two-qubit states of 
a given purity, the Werner state possesses the maximum 
discord. Next, we investigate the efficiency of dynamical 
decoupling schemes in preserving discord. 



B. Discord under Dynamical Decoupling 

Dynamical decoupling (DD) is a method of preserving 
coherences in NMR by frequent modulation of the spin- 
environment interaction with the help of a series of tt 
pulses 30 33] . The 7r-pulses effectively change the sign of 
the linear spin-bath interaction, and suitable time spac- 
ings between them undo the time evolution and suppress 
the T2-type decoherence due to the bath. We applied DD 
sequences immediately after obtaining the Werner state 
as in Eq. ([SS)) and followed that up with tomography. 
The CPMG DD involved a series of uniformly spaced tt 
pulses separated by 4-ms delays. For comparison, we la- 
bel as no-DD the evolution with the delays but without 
the pulses. The Uhrig DD (UDD) involved cycles of a 




FIG. 4: (a) Fidelity relative to the Werner state of the experi- 
mental and the simulated states as a function of 6, and (b) cor- 
responding discord values in units of e^/ln2. The maximum 
discord is obtained for the delay parameter 9 = (2n -I- l)vr/2, 
corresponding to preparation of Bell-diagonal states. 
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FIG. 5: (a) Fidelity of the experimental state relative to the 
Werner state for various DD schemes, and (b) corresponding 
discord values in units of e^/ In 2. 



seven pulse sequence [33|, |3J] . The time-instant tj of the 
j*^ TT pulse in each cycle was chosen according to Uhrig's 
formula tj — 28sin^(7rj/16) ms [33l |. 

Figure [5ja) shows the time dependence of fidelities rel- 
ative to the Werner state for no-DD, a CPMG DD se- 
quence, and UDD sequence. Figure [S^b) displays the 
corresponding discord values obtained using the exten- 
sive measurement method. We observe that the DD se- 
quences help in protecting fidelities of two-qubit quan- 
tum states, in agreement with an earlier work j34l |. On 
the other hand, there is not much difference between no- 
DD and DD schemes in preserving discord. We believe 
that the reason is the decay of purity during the DD 
sequences. This decay is a Ti-type decoherence, which 
the DD sequences do not protect against. By definition, 
fidelity looks at the orientation of the quantum state rel- 
ative to a target state and not its normalization. On the 
other hand, discord depends on the normalization of the 
state p. Our experiments indicate that though the DD 
sequences help prevent the quantum state from evolving 
to other quantum states, they are not useful in keeping 
the purity from decaying. 



C. Discord in Long-lived Singlet States 

Here we considered a pair of nuclear spins of the same 
isotope, i.e. Sample 2 shown in Fig. IHl^b). The closeness 
of their Larmor frequencies allows a spin-lock procedure 
[351] , where a continuous low-amplitude resonant rf irra- 
diation makes the spins nutate and drives the spin state 
toward an isotropic distribution with all Iz values equally 
likely. The pulse sequence for the preparation of a long- 
lived singlet state is shown in Fig. |6l[a). The rf wave for 
the spin-lock had the carrier frequency (uja + wb)/2, the 
mean Larmor frequency of the two spins, and a nutation 
frequency of 2 kHz. 

The average Hamiltonian in the interaction frame dur- 
ing the spin-lock interval is J^mt = 27rfi,J(/i ■/2), where /i 



and I2 are the spin operators. The singlet state, and the 
degenerate triplet states, form an orthonormal eigenbasis 

of this J^int'- 



\So 



1 



V2 
^i) = |00), iTo 



(|oi>-|io)), 

1 



V2 



(|oi) + |io)), |r_i 



(30) 
111). (31) 



As described earlier, the NMR system under ordinary 
conditions exists in a highly mixed state with a small 
purity. Leaving out the uniform distribution, the ground 
state is the Werner state, which is also called the long- 
lived singlet state (LLS) [3g, [S^. The LLS is antisym- 
metric with respect to spin exchange, and is not con- 
nected to other eigenstates (i.e. symmetric triplet states) 
by any symmetry preserving transformations such as the 
nonselective rf pulses and the intrapair dipolar interac- 
tion. Therefore, the LLS can survive for durations much 
longer than other non-equilibrium spin states. The LLS 
have been used in NMR experiments for the study of slow 
diffusion 13811 . for ultraprecise measurement of scalar in- 
teractions [39| , for storage and transport of parahydrogen 
[40l |4]| , and for preparation of high fidelity Bell states 
and other pseudopure states |43 |. 

The rf pulses prior to the spin-lock prepare a state 
which is a mixture of the \Sq) and |To) states, 

p{0) = ll + ii\So){So\-\To){To\). (32) 

During the spin-lock, rf irradiation mixes various com- 
ponents of states with the same spin, and the \To) state 
rapidly equilibrates with the other triplet states. Fur- 
thermore, all other coherences created due to pulse im- 
perfections also decay towards the background [ij]. On 
this equilibration, which takes a few seconds, the system 
reaches the Werner state. 



Plls 



l-(C/3), , e 



|^o)(^o|=PM/(e-e/3). (33) 



90x 180 




FIG. 6: The pulse sequence for preparing long-lived sin- 
glet states (a), molecular structure of 5-chlorothiophene-2- 
carbonitrile (b), and traceless real parts of the theoretical (c) 
and the experimental (d) density matrices. The experimental 
Werner state in (d) was obtained with a spin-lock of 16.4 s 
and has a fidelity of 0.99. 



This Werner state has a purity that differs from the one 
prepared from Sample 1, i.e., Eq. ([28| . 

To study the evolution of the density matrix state dur- 
ing spin-lock, we applied the spin-lock for a variable du- 
ration T and then carried out tomography to measure the 
traceless part of the density matrix p{t). We find that 
p{t) gradually evolves towards the Werner state pw, re- 
mains in that state for several tens of seconds, and ulti- 
mately decays towards the uniform state 1/4 that is the 
asymptotic eigenstate of the spin-lock evolution after the 
decay of all spin correlations. 

We monitored fidelity of the experimental state relative 
to the Werner state at 17 spin-lock durations, t = 2" ms 
with n = {0, 1, • • • , 16}). As shown in Fig. [71[a), it starts 
with a value of 0.85, reaches a maximum of 0.99 after 
a few seconds, and then decreases. The real parts of 
the deviation density matrices Jj^^ and p(t), with r = 
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uated fidelity, as shown in Fig. [Tja), we determined 
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FIG. 7: (a) Fidelity of the experimental state relative to the 
Werner state as a function of the spin-lock duration r, and 
(b) the corresponding values of discord in units of e^/ln2 
and geometric discord in units of e^/2. Discord values were 
obtained using the methods described in Sec. II. 



16.4 s corresponding to the niaximuni fidelity 0.99, are 
compared in Figs. Hfc) andHfd). 

Although fidelity is a good measure of how close a test 
density matrix is to the target density matrix, it does 
not capture the global decay of the purity of the density 
matrix. To monitor the decay of the purity as well as the 
closeness of the traceless parts of the density matrices, 
we use attenuated fidelity |25| . 



Fa^ 



Tr[j9(r) ■ p^] 
v/Tr[p(0)2]Tr[p2^] 



(34) 



It differs from fidelity in normalization, as evident from 
denominators of Eq. ([29l) and Eq. (|34| . and decreases as 
the purity p(t) decays. Figure [7l[a) also displays attenu- 
ated fidelity as a function of the spin-lock duration. We 
observe that it remains close to its initial value 0.85 until 
about 1 s, and then drops down. In particular, it starts 
dropping before the fidelity reaches its maximum value, 
and is 0.36 at r = 16.4 s. 

To understand the evolution of the density matrix 
during spin-lock, we consider a model consisting of ex- 
ponential equilibration of the |To) state with the other 
triplet states as well as an overall exponential decay of 
the singlet-triplet mixture toward the uniform identity 
state: 



1 



P{t) = 4I 



e"'^* I IPSO 



-Alt , 



-Ait\ 



Vto - (1 - e-^^^')pT^ 

(35) 

Here pr = \ m){T^\ + \T^){n\ + |r-i>(T_i|}, and Ai 
and A2 are the decay constants. By fitting the atten- 
uated fidelity of this model to the experimental atten- 



AT 



0.75 ms and A, 



26 s. These values indicate 



the rapid equilibration of the triplet states and the long- 
lived nature of the singlet state. It can be noticed that 
X^ is comparable to the rf period during the spin-lock, 
and X2 is significantly longer than spin-lattice relaxation 
time constants (Ti values). 

We also measured discord during the spin-lock evo- 
lution using the methods described in Sec. II, and 
the results are plotted versus the spin-lock duration in 
Fig. [Tljb). The results for the extensive measurement 
method D{B\A) and geometric discord DG{B\A) essen- 
tially agree (when scaled by appropriate factors), as ex- 
pected for accurate methods, and indicate that measure- 
ment errors in our experiments are rather small. By 
looking at angular variation of J{A : B) in the exten- 
sive measurement method, at r = 16.4 s when the state 
is closest to the isotropic Werner state, we estimate that 
the imperfections in our prepared LLS state give around 
3% error to discord values. 

The discord Dbd{B\A) obtained by assuming that the 
state is of Bell-diagonal form is an overestimate initially 
but becomes almost the same as the other two deter- 
minations beyond r = 1 ms. That means that artifact 
off-diagonal coherences are present in our prepared state, 
but they decay rapidly on a time scale comparable to 
Aj~^. The discord value for the two-parameter model of 
Eg. ([55]) is also shown in Fig.[71[b). It is accurate once the 
state becomes Bell-diagonal, but is unable to model the 
initial behavior. The reason for the initial discrepancy is 
that the off-diagonal components missing from Eg. ([55)) 
alter both I [A : B) and J {A : B). Later evolution and 
the asymptotic vanishing of discord after long durations 
of spin-lock is governed by the time scale A^ . 



IV. CONCLUSIONS 

Over the years, several hypotheses have been proposed 
to characterize the nonclassical part of quantum correla- 
tions, e.g., entanglement, violation of Bell-type inequal- 
ities, uncertainty relations, quantum discord, and so on. 
They are not all equivalent, especially for mixed states. 
Quantum discord aims to be the most inclusive of all of 
these, and we have focused our attention to the study of 
its dynamics in this work. 

It is well known that under ordinary NMR conditions, 
the purity of an ensemble of nuclear spin systems is below 
the threshold to exhibit entanglement. However, success- 
ful demonstrations of NMR quantum information pro- 
cessing indicate that quantum correlations do exist in 
such ensembles. In our work, we have, first, revisited 
the theoretical basis of discord as well as geometric dis- 
cord and then have described how they can be produced, 
manipulated, and monitored for an experimental density 
matrix, using the Werner state as a specific example. 
We have also implemented several checks in our investi- 
gations to keep track of experimental imperfections. 



The experimental study of discord was carried out in 
two different systems. In one system, we studied prepa- 
ration of discord using an entangling pulse sequence, 
and evolution of discord under dynamical decoupling se- 
quences. Under DD, discord did not show much improve- 
ment, but there was clear improvement in fidelity, which 
implies that DD protects against T2-type decoherences 
but not against Ti-type decoherences. In the second sys- 
tem, we used the construction of the long-lived singlet 
state to prepare the Werner state and examined its evo- 
lution under rf spin-lock. We can describe the accom- 
panying evolution of fidelity and discord reasonably well 
using a simple relaxation model. In both systems, the 



experimentally observed behavior and the fit parameters 
are in good agreement with expectations from theory and 
simulations. 
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